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Neutron resonance imaging (NRI) has recently emerged as an appealing technique for neutron radiography. 
Its complexity surpasses that of conventional transmission imaging, as it requires a high demand for both a 
neutron source and detector. Consequently, the progression of NRI technology has been sluggish since its in- 
ception in the 1980s, particularly considering the limited studies analyzing the neutron energy range above 
keV. The white neutron source (Back-n) at the China Spallation Neutron Source (CSNS) provides favorable 
beam conditions for the development of the NRI technique over a wide neutron energy range from eV to MeV. 
Neutron-sensitive microchannel plates (MCP) have emerged as a cutting-edge tool in the field of neutron detec- 
tion owing to their high temporal and spatial resolutions, high detection efficiency, and low noise. In this study, 
we report the development of a ‘°B-doped MCP detector, along with its associated electronics, data processing 
system, and NRI experiments at the Back-n. Individual heavy elements such as gold, silver, tungsten, and in- 
dium can be easily identified in the transmission images by their characteristic resonance peaks in the 1-100 eV 
energy range; the more difficult medium-weight elements such as iron, copper, and aluminum with resonance 
peaks in the 1-100 keV energy range can also be identified. In particular, results in the neutron energy range of 
dozens of keV (Aluminum) are reported here for the first time. 
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I. INTRODUCTION 


Neutron resonance imaging (NRI) is a nondestructive tech- 
nique for detecting the nuclide composition and distribution 
in samples [1—4]. Compared with traditional X-ray and neu- 
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clide composition and distribution in samples by using the 
characteristic resonances of neutron interactions. Imaging 
quality heavily relies on the neutron source and the perfor- 
mance of the neutron detector. Specifically, the detector is 
expected to achieve high temporal and spatial resolutions, a 
wide response range of the neutron energy, and a high detec- 
tion efficiency to obtain the energy-counting rate spectrum at 
each spatial pixel. A neutron source must have a wide energy 
spectrum, high flux, and good temporal resolution. There- 
fore, the NRI technique presents greater challenges, particu- 
larly for medium-energy neutrons ranging from hundreds of 
eV to hundreds of keV. Current neutron detectors have a low 
detection efficiency in this energy range, which complicates 
developing NRI to cover light to heavy elements. 


Many neutron-detector types, such as scintillator arrays, 
multi-pole photomultipliers, and position-sensitive gas detec- 
tors, are limited by their detection efficiency and spatial res- 
olution, thus cannot fully meet the requirements of NRI [10-— 
14]. CCD/CMOS cameras with gated time windows have 
been used in NRI to identify nuclides with resonance peaks 
in the eV energy range [15-17]. However, their low detec- 
tion efficiency and poor signal-to-noise ratio hinder attaining 
favorable results above the keV range. Neutron-sensitive mi- 
crochannel plates (nMCPs) have excellent properties includ- 
ing good temporal and spatial resolutions, irradiation resis- 
tance, a high gain, and low noise. They have been devel- 
oped for general neutron imaging, including thermal neu- 
tron imaging, and are considered suitable for the NRI tech- 
nique owing to their potential for a high neutron detection 
efficiency [18, 19]. The most commonly used nMCPs are 
Boron- and Cadmium-doped detectors. The nMCP matched 
with Timepix3 electronics has been used in NRI and demon- 
strated a good performance in the ultra-thermal neutron en- 
ergy region [20-24]. However, it is difficult to identify nu- 
clides with resonances above the keV range. As an analog- 
digital hybrid chip, the irradiation resistance of Timepix3 is 
affected by high-flux neutron sources. A suitable neutron 
source with a good time resolution and high flux covering a 
large energy range is rarely found. The white neutron source 
(Back-n) at the China Spallation Neutron Source (CSNS) was 
built in 2018 and serves as a multi-purpose platform to sup- 
port nuclear data measurements, detector calibrations, and 
irradiation-effect tests [25—28]. Back-n exhibits the advan- 
tages of a wide energy spectrum (0.5 eV-200 MeV), high 
flux (2x 10" n/cm?/s at 55 m from the target with P60 mm 
spot), low background, and good temporal resolution (<1%), 
which provide favorable beam conditions for exploiting the 
NRI technique over a wide neutron energy range from eV to 
MeV [29-32]. 


This study reports the development of a 1°B-doped MCP 
(called B-MCP hereafter) detector system that includes 
crossed anode strips, readout electronics, and data process- 
ing. Anode strips are suitable for NRI owing to their good 
performance in terms of the irradiation resistance, spatial res- 
olution (100 um), and counting rate (10 MHz/cm?). Several 
NRI demonstration experiments were conducted at Back-n to 
identify various elements with resonances in different energy 
regions, ranging from eV to dozens of keV, the performance 


of which exceeded that reported in literature. 


I. EXPERIMENTAL SETUP 


The detector system mainly consisted of a B-MCP/MCP 
stack, crossed anode strips, a front-amplifier module (FAM), 
a vacuum chamber, readout electronics (including digitization 
modules (FDMs) and a data aggregation module (DAM)), and 
a data acquisition computer, as shown in Fig. 1. 
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Fig. 1. Schematic of the B-MCP detector system. 


The B-MCP/MCP stack, as the core part of the detector, 
was composed of one B-MCP and two conventional MCPs 
encapsulated in a ceramic fixture. The B-MCP was devel- 
oped in cooperation with North Night Vision, and its main 
parameters are listed in Table 1. It contained approximately 
10 mol% !°B, which was owing to glass substrate doping, 
and the signal gain was approximately one order of magni- 
tude lower than that of the conventional MCP, which is bene- 
ficial for neutron detection. A large gain in the B-MCP leads 
to an excessive amplitude range of neutron signals, which is 
not conducive to the efficient acquisition of neutron events 
by specific electronics. The MCP pores had a bevel angle of 
approximately 8° and were stacked in a V-shape of two adja- 
cent MCPs to improve the secondary electron efficiency and 
reduce ion feedback. The voltage ratio of the B-MCP and the 
MCPs was 2:3:3 to further reduce the B-MCP gain, that is, 
the neutron signal amplitude range. 


Table 1. Main parameters of the B-MCP and a conventional MCP. 


Parameters B-MCP Conventional MCP 
Diameter(mm) 56 56 
Thickness(mm) 0.51 0.48 
Body resistance(MQ) 42 46 
Pore inner diameter(jzm) 8 10 
Pore area ratio (%) 60 60 
Gain @ 800V 120 2800 
10B (mol%) 10 = 


The a or “Li produced from the nuclear reaction of !°B(n, 
a)’Li with an incident neutron that can penetrate the sub- 
strate enters one or more pores of the B-MCP and excites 


the primary electrons on the inner wall. The electrons multi- 
ply many times in the pores to form an electron cloud at the 
exit surface of the MCP and were collected by the crossed 
anode strips to form the primary signal of the neutron case. 
The crossed strips consisted of diamond pads distributed on 
the surface layer of the PCB board. The pads were inter- 
leaved in both the X- and Y-directions, forming a sensitive 
area of 65x65 mm? with a total of 256 channels (the X- and 
Y-directions were both 128 channels). The strip pitch in both 
directions was 0.508 mm with a 0.10 mm gap between adja- 
cent pads. A schematic of this structure is shown in Fig. 2a. 
The gap between the exit surface of the MCP and the layer of 
the strips was 1.1 mm, allowing the electron cloud to "fire" 
multiple adjacent strips to facilitate the localization of the 
neutron case. 
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Fig. 2. (a) Local schematic structure of the crossed anode strips [33]; 
(b) physical object of the FAMs packed with the B-MCP/MCP stack. 


Primary signals were amplified and filtered by the FAMs 
before further processing, and the FAMs were packed with 
a B-MCP stack inside a vacuum chamber to avoid noise. 
The physical objects of the package are shown in Fig. 2b. 
The electronics system was developed in collaboration with 
the University of Science and Technology of China (USTC). 
Eight FDM boards outside the chamber were responsible for 
the waveform digitization of the amplified signals. The sam- 
pling frequency was set to 80 MHz to accommodate a signal 
pulse width of ~700 ns. The data were aggregated onto a 


DAM board for storage and offline analysis. The electron- 
ics system was triggered by a TO signal from a proton beam 
monitor (a fast current transformer FCT on the accelerator) as 
the starting time of the neutron flight, and the threshold was 
set to 50 fC to filter the noise of the electronics and wait for 
the signal trigger. The sampling duration was set to 0.6 u s 
before and 1.0 ws after the signal trigger, and the waveform 
on each strip was recorded. 
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Fig. 3. Experimental layout of the detector at the Back-n white neu- 
tron source. 


The experimental layout of the B-MCP detector at the 
Back-n white neutron source is shown in Fig. 3. Eight FDM 
boards were assembled in the chamber to optimize space- 
utilization and avoid direct beam irradiation. The surface of 
the B-MCP was approximately 76 m distant from the spalla- 
tion target. The beam spot at the detector was ®30 mm, which 
was smaller than the diameter of the B-MCP, allowing the 
neutron beam to completely pass through the sensitive area 
of the detector. The neutron flux was approximately 8.5 x 10° 
n/cm?/s, and the energy range was 0.5 eV-200 MeV [34]. 
The neutron pulse repetition rate was 25 Hz and had a double- 
bunch structure with an interval of 410 ns and a width of 60 
ns for each bunch. This introduced an error of more than 10% 
in the temporal resolution of neutrons above the MeV range, 
whereas for neutrons below dozens of keV, this effect was less 
than 1% and could be ignored in the experiments. 


HI. EXPERIMENT PROCESS AND DATA ANALYSIS 
METHOD 


An empty period of 6 ys after the TO trigger was set to 
avoid data blocking in the electronic channels owing to the 
complex responses of MCP components to the high-energy 
neutrons. Accordingly, the available energy ranged from 0.5 
eV to approximately 0.5 MeV. The amplitude spectra of the 
MCP dark counts and neutron signals were measured by in- 
creasing the voltage on the MCPs. Both the amplitudes of the 
dark count and neutron signal increased, and the two spectra 
gradually separated as the voltage increased. An appropriate 
voltage of 2400 V was selected to obtain a good neutron/dark 
count ratio, as shown in Fig. 4. Certain large neutron signals 


were saturated at an ADC channel of approximately 3500 ow- 
ing to the limitation of the electronics dynamic range. Con- 
versely, certain small signals were mixed in the dark counts 
and were difficult to separate. This is because the neutron sig- 
nals have a significantly wide amplitude range owing to the 
B-MCP gain. We set a second threshold at an ADC channel 
of approximately 900 to exclude the dark counts and noise. 
Certain small signals, whose amplitudes cannot be separated 
from the dark counts, will inevitably be lost. This reduces the 
overall detection efficiency, but does not affect the trend and 
details of the neutron spectrum because the amplitude of the 
neutron signal with an energy below a few hundred keV is 
nearly independent of the neutron energy. 
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Fig. 4. Amplitude spectra of the neutron cases and MCP dark counts 
at a voltage of 2400 V. 


Multiple adjacent channels are typically fired on the X and 
Y strips at each neutron event; a typical signal is shown in 
Fig. 5. Normally, the signal amplitude in the middle channel 
is higher than those on both sides, conforming to the elec- 
tron cloud distribution. A data analysis extracts the time-of- 
flight and hit positions of each neutron to plot the neutron 
energy spectrum and beam spot image. The stop time of the 
neutron flight is obtained using the peak of the unsaturated 
channels because the peak position has a stable amplitude; 
furthermore, unsaturated channels are always present, even 
for the saturated case. Peak deviations of more than two sam- 
pling points (approximately 25 ns) are considered different 
events within a sampling window, although these situations 
are significantly rare in the energy region below hundreds of 
keV. Both gravity and the center of the electron cloud dis- 
tribution can provide a spatial resolution of better than 100 
jum, which exceeds the width limit of the anode strips. How- 
ever, this significantly reduces the counts on each pixel and is 
not conducive to the neutron spectrum analysis for resonance 
imaging. Because the present efficiency of B-MCP is not suf- 
ficiently high, we used the central anode strip of the response 
strips in both the X and Y directions as the hitting point with 
a spatial resolution of 500 um depending on the pitch of the 
anode strips. 

Fig. 6a presents the neutron energy-count spectrum mea- 
sured on the Back-n source for 7200 s after deducting the 
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Fig. 5. Waveform of a typical neutron signal in the (a) X and (b) Y 
directions. 


dark counts and noise. The kinetic energy of neutrons can be 
calculated using Eq. (1). 


En = (y —1)mne? 


1 (1) 
Sa a 


where v = Lgtignt/TOF. The neutron flight length 
L flight 18 76 m from the neutron target to the surface of the 
B-MCP. The curve in Fig. 6a shows the response spectrum of 
the 1°B doped in the B-MCP at the Back-n source. Numer- 
ous detailed peaks ranging from 0.5 eV to tens of keV arise 
from the intrinsic microstructure of the Back-n energy spec- 
trum. The resonance peaks of certain nuclides such as 18!Ta 
(~4.2 eV), 9°Mo (~45 eV), °°Mn (~340 eV), and °°Fe (~28 
keV) in the spallation target are marked. Fig. 6b presents the 
Back-n beam spot image measured by the B-MCP detector. 
The half-height, full-width profile of the image was fitted to 
a diameter of approximately 30 mm, which coincided with 
the beamline setting. The outer profile of the image corre- 
sponds to the size of the B-MCP, that is 56 mm in diameter. 
The energy spectrum and beam spot of the neutron source can 
be simultaneously obtained in a single measurement using a 


B-MCP detector. Moreover, with a sufficiently high detec- 
tion efficiency (or total counts), a more detailed analysis can 
be performed on the energy spectra in different spot areas or 
beam spots at different energy intervals. Here, the neutron 
energy-count spectrum and beam spot of Back-n is demon- 
strated to provide a comparison of the data for the neutron 
resonance imaging experiments with the samples presented 
in Sec. IV. 
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Fig. 6. Neutron energy-count spectrum and (a) beam spot of the 
Back-n source (b) measured by the B-MCP detector. The identified 
image is distinctive by a yellow—blue scale. 


IV. NRI EXPERIMENTS 


The B-MCP detector can provide the position and energy 
information of each neutron case, such that the data can be 
flexibly processed in resonance imaging experiments to effi- 
ciently obtain more valuable information. In this section, the 
heavy elements including In, Au, Ag, and W with resonances 
in the eV energy region and the medium-mass elements in- 
cluding Cu, Fe, and Al with resonances in the keV energy 
region were selected to study the resonance imaging effect 
of the B-MCP detector on the Back-n source. For a sam- 


ple with an unknown composition, there are two approaches 
for determining the nuclides and their distribution in the sam- 
ple: 1) Comparison of the neutron energy-count spectra in 
the full image with and without a sample to determine the nu- 
clide components; the neutrons in the resonance region were 
then selected for imaging to obtain the spatial distribution of 
the corresponding nuclide. 2) Observing the sample distribu- 
tion using the total neutron transmission imaging; the neutron 
energy-count spectra at the sample locations were then ana- 
lyzed for the nuclide composition. A combination of these 
two methods allowed for a significantly adequate composi- 
tional analysis of the samples. In particular, for a small num- 
ber of nuclides in a sample, the latter is more effective. 


A. Nuclides with cross-section resonances in the 1 100 eV 
regions 


In the eV-energy region, four typical elements including 
Au, Ag, In, and W were selected, and their dominant cross- 
sectional resonances [35] with natural abundances are shown 
in Fig. 7. These resonances are distributed in the en- 
ergy range of 1-100 eV, in which the cross-sections of +97 Au 
and ''°In were approximately 30,000 barns, and the cross- 
sections of "*’Ag and "°’W were approximately 10000 barns, 
indicating a strong resonance absorption of neutrons in this 
region. In addition, the B-MCP detector has a relatively high 
detection efficiency in this region owing to the larger cross- 
section of 10B, which is approximately a few hundred barns. 
Therefore, the identification and imaging of these nuclides in 
this energy region are relatively easy. 
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Fig. 7. Cross-sections of 115Ty ratw 197 Au, and nat Ag within the 
energy region of 1-200 eV. 


Fig. 8a presents the objects of the four samples in different 
shapes with the same thickness of 1.0 mm. These samples 
were placed at the center of the neutron beam, 70 mm distant 
from the B-MCP detector, such that the neutrons vertically 
passed through the sample plane. Using the aforementioned 
experimental conditions, the measurement was performed for 
3 h and the total-transmission imaging of the samples was 
obtained, as shown in Fig. 8b. Initially, we compared the 
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Fig. 8. (a) Sample objects and (b) the neutron transmission image over the full-energy region. 


total neutron energy-count spectra with and without these 
samples. The results show no apparent difference between 
these two spectra, and the strong resonance absorptions are 
invisible (the two identical spectra are the same as those in 
Fig. 6a, implying that the nuclide composition cannot be de- 
termined first. The number of these samples was relatively 
small, and the large amount of data in the sample-free re- 
gion served as a background to extract the resonance peaks of 
the samples. However, in the total-transmission images, the 
profiles of the four samples were clearly visible. The image 
underwent background subtraction and was normalized to the 
comparison data without the samples. The brightness reflects 
the intensity of neutron absorption, that is, the amount of the 
sample. It is impossible to identify the sample components 
in the total-transmission image without their shapes because 
the brightness value is related to both the cross-sections of the 
nuclides and their areal density on the pixels. 


We selected the neutron energy-count spectra from the 
sample areas and normalized them to the comparison data 
to obtain the ratio spectra with and without the samples, as 
shown in the ratio spectra in Fig. 9. The peak clusters in the 
ratio spectra reflect the nuclide species, and the intensity of 
the peaks reflects the areal density of the nuclides, which suf- 
ficiently corresponds with the primary and secondary peaks 
of the four nuclides, namely '!°In, 197 Au, "**W, and ”% Ag. 
These peaks act as fingerprints to identify the presence of 
nuclides in the samples. We selected the energy regions 
with these characteristic peaks for full-image neutron reso- 
nance imaging. As shown in Fig. 9, the four samples of 
115Tn, 197 Au, "@*W, and ”% Ag can be independently identi- 
fied using the spatial distribution. Each image was normal- 
ized and divided by a background image that selected the 
non-resonance-energy region adjacent to the dominant peak 
of the nuclide to represent each nuclide separately. This im- 
age processing was discussed in a previous article [17]. ''°In 
and "“*W can be easily identified owing to their distinct and 
non-overlapping resonances at approximately 1.4 eV-1.5 eV 
and 18.0 eV-19.5 eV, respectively. However, the dominant 


peaks of the 197 Au and ”% Ag overlap at approximately 4.5 
eV-5.5 eV, which causes them to simultaneously appear dur- 
ing resonance-identification imaging. To distinguish between 
197 Au and ”% Ag, we selected a part of the 1°97 Au peak at 
approximately 4.8 eV-5.0 eV to exclude "Ag, as well as 
the secondary peaks of "“* Ag at approximately 16.0 eV-16.5 
eV. This reduces the image statistics, which can be solved by 
using longer measurement durations. In the absence of the 
coexistence of 197 Au and ”‘Ag, the imaging effect using the 
dominant peaks was better. The nuclides with cross-sectional 
resonances in the eV energy regions can be easily identified 
and imaged using the B-MCP detector at the Back-n source 
with the advantage of strong resonance peaks. 


B. Nuclides with cross-section resonances in the 1-100 keV 
regions 


Fig. 10 presents the cross-sections of "“’Cu, "*’Fe, and 
27 Al over a wide energy range from 100 eV to 500 keV [35]. 
The resonance peaks were mostly located in this region, and 
the peak heights were nearly two orders of magnitude lower 
than those of the heavy elements in Sec. IV A, which compli- 
cates nuclide resonance imaging owing to a reduction in the 
imaging contrast. Moreover, the peaks were significantly nar- 
row, which further reduces the cross-sectional integral. Be- 
cause the widths of the resonance peaks in this energy re- 
gion were comparable to the spacing between the two neu- 
tron bunches, the pulse structure introduces approximately 
14%-50% of the undesired neutrons into the selected reso- 
nance ranges [17]. This effect can be resolved with sufficient 
statistics by using energy-resolved processing, in which the 
double-bunch structure can still be used to identify these ele- 
ments by selecting energy-resolved neutrons [36]. 

Fig. lla shows the Cu, Fe and Al samples with the same 
thickness of 10 mm. Fig. 11b shows the neutron transmis- 
sion image in the full-energy region, in which all the con- 
stituents are clearly visible. The Cu and Fe samples exhibited 
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Fig. 9. Identification imaging with the selected energy region cor- 
responding to the resonance peak of (a) 1151n, (b) +97 Au, (c) "tW, 
and (d)"*" Ag. 


high and similar brightness values owing to their strong neu- 
tron absorption. The Al area is the darkest because its cross- 
sectional integral over the full-energy region (including the 
nucleus number density) is smaller than those of Fe and Cu, 
as illustrated in Fig. 10. As expected, the individual nuclides 
were unable to be identified via total-transmission imaging. 
Resonance imaging analyses were performed to identify 
the samples and their spatial distribution by using the same 
method described in Sec. IV A. The ratio spectra with and 
without the samples in the sample areas are shown in Fig. 12. 
The peak clusters in the ratio spectra reflect the nuclide 
species, which sufficiently correspond with the resonances of 
the sample nuclides "°’Cu, "*’Fe, and 2” Al. We selected the 
energy regions with these characteristic peaks for full-area 
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Fig. 10. Neutron cross-sections of "*'Cu, "*’Fe, and ?’Al in the 
energy region of 100 eV—500 keV. 
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Fig. 11. (a) Sample objects and (b) the neutron transmission image 
over the full-energy region. 


resonance imaging, as shown in Fig. 12. The three samples of 
Cu, Fe, and Al could be independently identified with a spa- 
tial distribution, although not as well as the samples in the eV 
energy region. This is partly because the intensity and width 
of these resonances are not large and because the detection 
efficiency of '°B in the keV region is reduced. In particu- 
lar, for the Al sample, the three wider resonances at approxi- 


mately 35, 90, and 150 keV were selected instead of the high- 
est resonance. They have a cross-section of only 20 barns, 
but the signal-to-noise ratio is significantly better than that of 
a higher but narrow peak. The nuclides with cross-sectional 
resonances in the keV energy regions can be identified and 
imaged at the Back-n source with a good signal-to-noise ratio 
of the B-MCP detector. One advantage of the B-MCP detec- 
tor over other methods such as the CMOS/CCD camera is that 
it can simultaneously obtain multiple resonances of nuclides 
for superimposed imaging to improve the signal-to-noise ra- 
tio and resonant image contrast. 
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Fig. 12. Identification imaging with the selected energy region cor- 
responding to the resonance peaks of (a) "*’Cu, (b) "“*Fe, and (c) 
27 

Al. 


Ideally, B-MCP can provide neutron energy-count spectra 
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